Warm turbulence in the Boltzmann equation 
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We study the single-particle distributions of three-dimensional hard sphere gas described by the 
Boltzmann equation. We focus on the steady homogeneous isotropic solutions in thermodynamically 
open conditions, i.e. in the presence of forcing and dissipation. We observe nonequilibrium steady 
state solution characterized by a warm turbulence, that is an energy and particle cascade superim- 
posed on the Maxwell-Boltzmann distribution. We use a dimensional analysis approach to relate the 
thermodynamic quantities of the steady state with the characteristics of the forcing and dissipation 
terms. In particular, we present an analytical prediction for the temperature of the system which 
we show to be dependent only on the forcing and dissipative scales. Numerical simulations of the 
Boltzmann equation support our analytical predictions. 
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Three-dimensional fluid turbulence is characterized by 
an energy cascade process whose concept was introduced 
by Richardson in 1922 and developed by Kolmogorov in 
1941 [II [5]. The idea is that if one injects energy at some 
scale, creating eddies, then those eddies may interact and 
create smaller and smaller eddies up to the dissipation 
scale. The result is the famous power-law Kolmogorov 
spectrum which corresponds to a constant flux of energy 
from large to small scales. 

Some years later, in 1965 [3], Zakharov found that 
dispersive weakly nonlinear wave systems, e.g. surface 
gravity waves, may exhibit a very similar phenomenol- 
ogy as fluid turbulence, i.e. a constant flux of energy to- 
wards small scales, a direct cascade. Besides ocean waves 
[1], these cascades have been studied in a large num- 
ber of weakly nonlinear systems such as internal waves 
[5], nonlinear optics Bose-Einstein condensation [3- 
IH], magnetohydrodynamics [TU]. Such regime, valid for 
weak nonlinearity, is now known as weak wave turbu- 
lence and the power-law flux carrying states are called 
Kolmogorov-Zakharov (KZ) spectra. The extraordinary 
fact is that the KZ solutions corresponding to the tur- 
bulent cascades are found to be exact analytical solu- 
tions of a wave kinetic equation which describes the evo- 
lution of the wave spectrum. The wave kinetic equation 
has a structure that resembles the classical Boltzmann 
equation: the evolution of the distribution is driven by a 
collision integral which conserves mass, momentum and 
energy. 

Keeping in mind such an analogy, our aim is to address 
in the present Letter the following fundamental ques- 
tions. Let us consider a gas whose single-particle distri- 



bution evolution is described by the Boltzmann equation 
and suppose that we inject in the system particles with 
a specific energy and remove those that reach an energy 
larger/smaller than a preselected threshold. How does 
the distribution function change with respect to the ther- 
modynamic equilibrium solution? How the temperature 
and chemical potential will be affected by the presence 
of forcing and dissipation? 

In order to answer we will first consider numerical sim- 
ulations of the Boltzmann equation in the homogeneous 
and isotropic case and then present an argument based 
on dimensional analysis that allows us to explain the nu- 
merical results and make some predictions on the steady 
nonequilibrium properties of the system. Our starting 
point is the homogeneous, forced and damped Boltzmann 
equation 

-Q^ = I coll + F - D, with 

^+oo ^ ^ (1) 

/co/i = / Wi2^ [n[n2- nin2]dv2dv[dV2, 

where n,; = n(x,Vi,t) is the single-particle distribution 
function and primes denote particles after the collision. 
We have included a source term F and a sink D which 
will be specified once the numerics are discussed. As we 
consider elastic collisions, the general way to express W 
is 

W^^'"' =<jSivi+V2-v[~v'.,)S{Ei+E2-E[~E!,), (2) 

with Ei — |v,;p/2 being the kinetic energy per unit mass. 
The (5-functions assure conservation of the total momen- 
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turn and the total kinetic energy. In this Letter, we con- 
sider the three-dimensional rigid sphere gas {a is inde- 
pendent of v) in isotropic conditions. 

In absence of forcing and dissipation any initial condi- 
tion will relax to the Maxwcll-Boltzmann (MB) distribu- 
tion 
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nMB{E) = noe 



(3) 



where A = uq e^T with /i and T the chemical poten- 
tial and temperature of the system respectively. This 
equipartition mechanism, consequence of the H-theorem, 
has been checked as a benchmark of our numerical code. 
In order to consider an open system we have then in- 
cluded forcing and dissipation written in i?-space. The 
forcing term is constant in time and has the role of inject- 
ing particles with energies narrowly concentrated around 
some value Ef, i.e. F{E) = F if \E — E-f \ < 5f and zero 
otherwise, where _F is a positive constant. The dissipa- 
tion term is implemented as a filter which removes, at 
each iteration time, energy and particles outside of the 
domain (i^min, £^max), i-e. D{E) = if G (iJmin 
and D{E) = — oo otherwise. In such conditions, a di- 
mensional analysis shows that ([I]) may have also solutions 
characterized by a constant flux of energy and mass which 
corresponds to the KZ solutions. Indeed, the particle flux 
T] and the energy flux e can be estimated as: 



^ ^ 2^ / ^ E^/^dE ^ n^E^^ 

e = 2TT r —E^^^dE^n^E^/\ 
dt 



(4) 



Assuming that one of the fluxes is constant trough energy 
scales, we can immediately derive the KZ solutions ~ 
E~y* and n.'^E-^/'^ 31]. 

In our simulations the computational domain is uni- 
formly discretized in 501 points having AE = 1 and the 
initial distribution is n{E, t — 0) — 0. We flrst consider 
the case characterized by E'min = 5, E'max = 250 and 
the forcing located at energies between 35 and 37. Nu- 
merical results for the final steady states evaluated for 
three different forcing amplitudes, F = 10^"*, 10^^, 10^^, 
are presented in Fig. [T] Independently of the forcing, 
a stationary distribution is reached in the simulations, 
implying the existence of a flux from the forcing region 
to the boundaries of the domain. The steady solutions 
of our simulation do not appear to be very far from the 
Maxwell-Boltzmann distribution ([3|, which corresponds 
to a straight line in a lin-log coordinates. It is clear 
from the plot that the only consequence of increasing the 
forcing rate is to shift upwards the curves, leaving un- 
changed the slopes. This is the first indication that the 
temperature of the system remains constant, indepen- 
dently of the forcing. We can also compute, as a func- 
tion of time, the average energy per particle pe/ Pm, with 
PM = 27r / n{E) E^''^ dE and = 27r / n{E) E^'"^ dE. 



10 



■3 J 



10"' 



10" 




250 



FIG. 1: Final steady states of the Boltzmann equation ([T]) 
obtained for different values of the forcing amplitude F. Axes 
are in lin-log coordinates. -Emin = 5, iJmax = 250 and Ef = 
35 - 37. 
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FIG. 2: Time evolution of the average energy per particle for 
different forcing amplitudes. Axes are in log-lin coordinates. 

-Emin = 5, -Bniax = 250 and = 35 - 37. 



As shown in Fig. [2j for large times this quantity reaches 
a unique value for the three simulations considered. We 
recall that for a pure MB distribution such ratio is pro- 
portional to 3/2 the temperature of the system. 

We have not observed KZ solutions in our simulations. 
The reason is that the interactions are non-local in scales, 
as already pointed out in [T^] : the collision integral does 
not converge for such solutions. Moreover, the energy 
and particle flux directions, [13] , associated to such solu- 
tions, are opposite to the one predicted by the Fj0rtoft 
argument which imposes that the energy should have a 
direct cascade, i.e. from low to high energies, while par- 
ticle an inverse one [l4] . 

As shown in Fig. [l] we observe distributions not far 
from the thermodynamic equilibrium Q solution. How- 
ever, we are in a forced and dissipated situation and let us 
assume the existence of a small but flnite flux correction 
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"living on top" of the Maxwell-Boltzmann distribution. 
This behavior, named warm cascade, has aheady been 
observed in other physical systems O [15] and is charac- 
terized by constant flux cascades perturbing the thermo- 
dynamic equilibrium distribution. Mathematically, we 
assume that 

= ?^AfB (1 + ri), (5) 

with n the deviations with respect to the MB distribu- 
tion which are responsible for the fluxes; note that not 
necessarily n is small with respect to one. If we plug such 
ansatz in the equations for the fluxes Q we obtain 

V = ci nljB h{2 + n) E'l^ 

(6) 

where ci and C2 are two constants which cannot be de- 
termined through dimensional analysis. The term nj^jg, 
corresponding to the unperturbed MB distribution, does 
not give any contribution because it is not responsible for 
any net flux. 

Our aim is to relate the macroscopic properties of the 
system with the forcing and dissipation rates. From our 
numerical computation we observe that, as we get closer 
to the cut-off scales (i?min and i?max), deviations from 
a pure MB distribution becomes more relevant. Conse- 
quently, h becomes of the order one for E = -Emin or 
E = E'max and therefore from ([6| we have 

e = C2 A^e-'-^El(l, 

where we have redefined the constant ci and C2 . We now 
verify the above relations through the direct computation 
of the Boltzmann equation. In Fig. |3]we show the depen- 
dence of A on the fluxes for three simulations previously 
described (Fig. [T]). Supposing to use as temperature of 
the system the simple relation coming from the MB dis- 
tribution T = 2pe/5pm, we can observe in Fig. |3] (and 
its inset) that A scales as the square root of the incoming 
fluxes for fixed dissipative scales i?min and i^max- From 
the numerics we can estimate, by a fit, the constants Ci 
and C2 whose values are reported in the figure. 

We are now able to predict the dependence of the tem- 
perature on the forcing and the dissipative scales. As- 
suming that they are widely separated, that is E'min ^ 
Ef i^maxj we have e — TjEf. Then from equations ([T]) 
we get 

rp _ 2 (jj/niax ~ £'min) 

I In £:,„ax - i In -Emin -InEf + C3 

with C3 = In (c2/ci). The temperature of the system does 
not depend on the incoming fluxes but only on the forcing 
and dissipative scales; this is consistent with results of the 
numerical simulations presented in Fig. [l] and [2] 
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FIG. 3: ^ as a function of the fluxes. Circles (and squares 
in the inset) correspond to the numerical simulation results. 
The continuous and dashed lines are the predictions ([7| where 
the temperature, T — 2pE/3pM, has been defined using the 
pure MB distribution. The values of the constant Ci and C2, 
shown in the plot, are found by a fit. 

We can check the validity of prediction ([s]) by con- 
sidering different simulations of the isotropic Boltzmann 
equation ([T]) changing the forcing and dissipative scales. 
Before entering in the details, we emphasize that our pre- 
dictions are based on a dimensional argument and Ci and 
C2, which define C3, can only be measured via numerical 
computations. In our simulations we have observed that, 
as Emin, Ef and i?max changes, ci and C2 assumes dif- 
ferent values. This may be related to the fact that we 
are considering the hypothesis that the particle fiux is 
all dissipated at low energy scales and all energy flux 
at high ones, that is Emin <C -E/ ^ Emax- This is not 
always verified in the numerical simulations due to finite- 
ness of the computational domain. In our numerical sim- 
ulations we have measured an upper and lower bound 
for our constants leading to —4.21 < C3 < —1.05. This 
interval includes the analytical prediction C3 — 2 In (2/9) 
obtained using a diffusion approximation model in the 
limit Emin ^ T ^ E^ax presented in |14j . 

In Fig. ID we show the comparison between the esti- 
mation of the temperature from Q and the numerical 
simulation varying the low energy dissipation scale Emin. 
The incoming flux, the forcing and the high energy dissi- 
pative scale have been kept fixed to the respective values 
of E 10~^, Ef = 36, Sf = 1, Emax = 250. The dashed 
and the solid line in the figure are the predictions ob- 
tained with C3 = —1.05 and C3 = —4.21, respectively. 
The agreement between the temperature prediction and 
the computation is satisfactory. 

In Fig. [5] we present the behavior of the temperature 
as a function of the high energy cut-off scale Emax, with 
F = 10-^, E„iin = 5 and Ef = 36, 5/ = 1. The compari- 
son between our estimation ([s]) and the numerics is very 
good. 
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FIG. 4: Temperature as a function of the low energy cut- 
off Bmin, keeping fixed F = 10~^, Ej = 36, 5/ = 1 and 
£^max = 250. Circles represent the results from numerical 
computations. The continuous and dashed lines are the pre- 
diction Q for different values of C3. 
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FIG. 5: Temperature as a function of the high energy cut-off 
iSmax, keeping fixed F = 10~^, iSmin = 5 and Ej — 36, 5/ — 1. 
Circles represent the results from numerical computations. 
The continuous and dashed lines are the prediction Q for 
different values of C3 . 

In the present Letter we have investigated the station- 
ary states of a three-dimensional hard sphere gas whose 
single-particle distribution function is modeled by the ho- 
mogeneous isotropic Boltzmann equation. In particular, 
we were interested in the steady nonequilibrium states 
in an open system condition, i.e. in the presence of a 
forcing and dissipation mechanisms. Using the language 
of wave turbulence theory we have assumed the presence 
of a warm cascades as steady distributions. These solu- 
tions are characterized by constant particle and energy 
fluxes superimposed on the thermodynamic Maxwell- 
Boltzmann distribution. Our assumption has allowed us 
to relate the thermodynamic quantities of the system to 



the characteristics of the forcing and dissipative terms of 
the equation. 

In particular, we have been able to give an analytical 
relation for the temperature reached in the system as a 
function of the forcing and dissipative scales. One of the 
main results is that the temperature is independent of 
how much energy and particles we inject in the system 
but depends only on the cut-off scales and on the forcing 
scale. By a direct numerical integration of the Boltzmann 
equation we have shown that our numerical results are 
consistent with the aforementioned prediction. 

We believe that the approach used in this Letter may 
open some perspectives towards understanding nonequi- 
librium steady states in other physical systems via analo- 
gies to cases widely studied by the wave turbulence the- 
ory. 
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